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Abstract 

Background: Recent releases have been carried out with Aedes oegypti mosquitoes infected with the wMelPop 
mosquito cell-line adapted (wMelPop-CLA) strain of Wolbachia. This infection introduced from Drosophila provides 
strong blockage of dengue and other arboviruses but also has large fitness costs in laboratory tests. The releases 
were used to evaluate the fitness of released infected mosquitoes, and (following termination of releases) to test for 
any effects of wMelPop-CLA on wing size and shape when mosquitoes were reared under field conditions. 

Methods: We monitored gravid females via double sticky traps to assess the reproductive success of 
i/vMelPop-CLA-infected females and also sampled the overall mosquito population post-release using Biogent 
Sentinel traps. Morphometric analyses were used to evaluate infection effects on wing shape as well as size. 

Results: Oviposition success as assessed through double sticky traps was unrelated to size of released mosquitoes. 
However, released mosquitoes with lower wing loading were more successful. Furthermore, i/i/MelPop-CLA-infected 
mosquitoes had 38.3% of the oviposition success of uninfected mosquitoes based on the predicted infection 
frequency after release. Environmental conditions affected wing shape and particularly size across time in 
uninfected mosquitoes, but not in naturally-reared i/\/MelPop-CLA-infected mosquitoes. Although the overall size 
and shape do not differ between naturally-reared i/\/MelPop-CLA-infected and uninfected mosquitoes, the infected 
mosquitoes tended to have smaller wings than uninfected mosquitoes during the cooler November in comparison 
to December. 

Conclusion: These results confirm the lower fitness of i/i/MelPop-CLA infection under field conditions, helping to 
explain challenges associated with a successful invasion by this strain. In the long run, invasion may depend on 
releasing strains carrying insecticide resistance or egg desiccation resistance, combined with an active pre-release 
population suppression program. 
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Background 

The virulent maternally-inherited endosymbiotic Wolba- 
chia strain, wMelPop-CLA [1], reduces lifespan [2,3] and 
strongly blocks dengue virus proliferation in the primary 
dengue mosquito vector, Aedes aegypti [4-6]. Replacing 
field mosquitoes with wMelPop-CLA-infected individuals 
could therefore reduce dengue transmission in the field. 
However, wMelPop-CLA has deleterious effects on its 
host including reduced viability of mosquitoes in quiescent 
eggs [3,7], reduced fecundity [3,8], reduced ability to 
blood-feed [9,10], and altered development [3,11]. This 
poses significant challenges for the spread of the wMel- 
Pop-CLA infection in the field [3,12], unless such deleteri- 
ous effects can be overcome through combined modalities 
[13]. However, these deleterious effects by themselves 
might provide an effective method for population suppres- 
sion [3,7] particularly in genetically isolated sites [14]. 

Following the successful invasion of mosquitoes with 
the less virulent strain of Wolbachia, wMel in two isolated 
suburbs of Cairns, Northern Queensland, Australia in 
2011 [15], field trials releasing wMelPop-CLA in two other 
regions were undertaken in 2012 [16]. These releases pro- 
vided an opportunity to assess the success of wMelPop- 
CLA-infected females in sourcing an oviposition site and 
to assess morphometric traits to evaluate fitness of the re- 
leased mosquitoes. This follows earlier work [11] showing 
that larger females tend to be more successful at finding 
oviposition sites, and suggesting that infected mosquitoes 
emerging in the field have similar sizes to uninfected mos- 
quitoes, at least for the non-virulent wMel strain. 

We monitored fitness based on successful oviposition 
and morphology of wMelPop-CLA-infected mosquitoes 
during the first three weeks of the release, and also geo- 
metric morphology post-release at the transition from 
the dry to wet season. The following questions were 
considered. (1) Are there differences in the measure- 
ments of morphometric traits between field-destined 
and recaptured released wMelPop-CLA-infected females? 
(2) Is the frequency of successful ovipositing females 
that are infected and released comparable to the estimated 
release frequency? (3) Will naturally-reared infected 
mosquitoes have differing body size and shape in con- 
trast to uninfected mosquitoes? We then use these find- 
ings to help interpret the outcome of the wMelPop-CLA 
releases in terms of whether the released females are 
surviving, mating, blood feeding and ovipositing, as well 
as the fitness of subsequent generations in the field. For 
this purpose, we employed two mosquito trapping sys- 
tems: double sticky traps (DST) [17] to monitor oviposit- 
ing female mosquitoes at the beginning of the release 
and BioGent Sentinel traps (BGS) to monitor the adult 
mosquito population several months after release [18-20]. 
We obtained body size and shape measurements from 
mosquitoes caught in both types of trap. 



Methods 

Released mosquitoes 

The mosquitoes came from releases undertaken in 2012 at 
Machans Beach near Cairns in northern Queensland [16]. 
Releases were undertaken with offspring of wMelPop- 
CLA-infected mosquitoes established in two semi-field 
cages (8.0 m x 9.0 m x 4.1 m) [21] following a similar pro- 
cedure described in Hoffmann et al [15] where uninfected 
males derived from the field are periodically released into 
the cage to counter the effects of inbreeding. Offspring are 
established from eggs laid by these females and released 
soon after emergence. Releases took place over a 13-week 
period at Machans Beach and involved around 10 females 
being released weekly per house. Fitness assessments for 
released individuals were based on the first 10-day period 
of the release when individuals from multiple releases 
would not have been captured in traps (see below). 
Although these releases ultimately did not result in the 
wMelPop-CLA infection becoming established in Machans 
Beach [16], the infection persisted for several months 
after releases were terminated and this provided an op- 
portunity to investigate fitness of the release material as 
well as morphometric effects of the infection in mosqui- 
toes reared under field conditions. 

Morphometric traits and fitness of released mosquitoes 

One hundred DSTs were employed in fifty residential lo- 
cations (two traps each) in Machans Beach in the first 
week of the mosquito release. The DSTs [17] were made 
up of two black containers: the bottom container held 
1L of water and a lucerne pellet infusion [22], while the 
top was clipped to the bottom and contained a sticky 
panel (UVR-32, Atlantic Paste and Glue) on the internal 
surface and a gap on top to allow mosquito entry. As a 
DST contains an infusion attractive to Aedes, up to 99% 
of the trapped A. aegypti consist of gravid/parous fe- 
males [23]. Virgin gravid females are less receptive to 
oviposition site seeking [24,25], and females caught in 
DSTs were, therefore, assumed to be mostly mated. 
Samples were collected on 4th, 6th, 9th, 11th, 13th and 
16th January 2012 and tested for Wolbachia infection. 

Measurements of wing centroid size, thorax length 
and the wing size to thorax length (wing/thorax) ratio 
were compared among samples involving field-caught 
released mosquitoes from the DSTs, infected mosquitoes 
from the field cage and uninfected mosquitoes from the 
DSTs. Coefficients of variations (CoV) associated with 
the three morphometric measurements and wing shape 
(see below) were also compared. 

Comparing ovipositing infected females with estimates 
based on release frequencies 

As releases of infected mosquitoes began on 4th January 
2012 and these mosquitoes had not been blood-fed, 
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wMelPop-CLA-infected females would not have been 
gravid on 4th January and 6th January. Females require 
approximately 2 to 3 days at 26-28°C to develop mature 
eggs post-blood feeding prior to oviposition [24,26,27]. 
Thus, only infected females caught after these dates were 
considered. 

Assuming an equal and constant rate of trapping of 
non-gravid females for both uninfected and wMelPop- 
CLA-infected females, we estimated the infection fre- 
quency of ovipositing females directly attributed to the 
first release. We estimated the proportion of ovipositing 
infected females in the first release using data from col- 
lections on 9th, 11th and 13th January, prior to females 
from the second release contributing to traps. 

We also estimated the expected number of mosquitoes 
in their first gonotrophic cycle from the first release to 
be caught in the traps after 13th January. This provided 
an estimate of the proportion of successfully ovipositing 
females within the first 10 days of release (4th- 13th 
January). We did this by fitting two models. In the first, 
we fitted an exponential decay function on the daily 
trapping rate of wMelPop-CLA-infected females with re- 
spect to time. This was done via linear regression of the 
natural log of the daily trapping rate against the mid- 
point date between collection dates, excluding the first 
two, and the last collections. The second model was a 
logistic growth model on the cumulative number of 
mosquitoes trapped over time. 

We assumed negligible contribution from the second 
gonotrophic cycle due to reduced survival of the wMel- 
Pop-CLA-infected mosquitoes [2]. Released wMel-in- 
fected females are likely to have a field survival of 70-90% 
per day [28] and around 5-30% are expected to survive to 
the second gonotrophic cycle (10-15 days old). The sur- 
vival rate is expected to be somewhat lower for wMelPop- 
CLA-infected females. 

Effects of infection on size and shape in generations 
subsequent to release 

To examine the effect of the infection on traits in 
naturally-reared mosquitoes, samples from November- 
December 2012 were collected to study size and shape 
of mosquitoes at a time when the population remained 
polymorphic for the wMelPop-CLA infection, but well 
after releases had been terminated. BGS -traps were 
placed at around 100 residential properties with the con- 
sent of the owners. Traps were inspected once a week. 
BGS -traps are effective at capturing all adult stages of A. 
aegypti mosquitoes including young nullipars [18-20]. 

Wing size and shape were compared in mosquitoes 
from the BGS -traps; it was not possible to obtain thorax 
measurements because these were used immediately to 
screen for Wolbachia infection. Comparisons were made 
for wing size, CoV and shape between uninfected and 



infected groups, and also between collections. We used 
the total number of mosquitoes caught in the traps as 
an indicator of mosquito density and also noted the 
average temperature during the 20 days leading up to 
trap collection. 

l/l/o/bach/a-infection status 

The DST samples were tested using the protocol detailed 
in Lee et al [29] (see also Additional file 1: Table SI). PCR 
conditions proceeded with the following settings: 95°C for 
10 minutes, 40 cycles of 95°C for 5 seconds, 58°C for 15 
seconds and 72°C for 15 seconds, ending with a 95°C 
1 -minute heating followed by cool down to 40°C for 20 
seconds before raising to 65°C. A melting curve analysis 
was performed immediately via a gradual increase of 
temperature from 65°C to 95°C. Using a Light Cycler 480 
(Roche Applied Science), we determined the crossing 
point (Cp) values and melting temperatures (Tm) which 
allows us to determine presence or absence of Wolbachia. 
The Wolbachia specific primers in Lee et al were replaced 
with wMelPop-CLA specific primers initially developed to 
screen the BGS-trap samples (see below). 

The BGS-trap samples were screened using a different 
protocol. Mosquitoes collected in BGS traps in North 
Queensland were kept in 70% ethanol and shipped to 
Monash University, Melbourne for PCR. For DNA extrac- 
tions, adult mosquitoes were washed in milli Q (Millipore) 
water and individually transferred to 96-well PCR plates 
containing a 2 mm glass bead and filled with 50 uL extrac- 
tion buffer per well. The extraction buffer consisted of 4.8 
mL squash buffer (10 mM Tris pH 8.2, 1 mM EDTA, 50 
mM NaCl) and 60 uL proteinase K (15 mg/ml, QIAGEN). 
Each extraction plate included wMel, wMelPop-CLA and 
uninfected mosquitoes from laboratory colonies as controls. 
After homogenizing the samples in a Mini-Beadbeater 
(Biospec Products) for 1.5 min, the plates were incubated 
in a 96-well thermocycler block for 5 min at 56°C then 
boiled at 95°C for 5 min (to inactivate proteinase K) and 
cooled down to 4°C. Extracted DNA was then stored at 4°C 
for up to 3 days prior to qPCR. 

The primers and probes are described in Additional 
file 1: Table SI. The A. aegypti primers and probe se- 
quences were designed for the rpsYl gene [4]. The 
primers and probes specific for wMel and wMelPop- 
CLA strains were designed using the PrimeTime qPCR 
tool (Integrated DNA Technologies). The wMel specific 
primers and probe were designed to target the WD0513 
gene that is present in wMel [30] but is absent in 
the wMelPop-CLA genome that was transinfected into 
the PGYP1 mosquitoes (Woolfit, unpublished results). 
The primers and probe specific for wMelPop-CLA were 
designed across the region spanning the IS5 element 
inserted into WD 13 10 [31] (and Riegler et al, unpub- 
lished data), since this ISS insertion is absent in the 
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wMe\ WD 13 10 gene [30]. Each of the three probes was 
labelled with a different fluorophore with non-overlapping 
emission wavelengths in order to be used in the same 
multiplex reaction: rpsYl with 6-carboxyfluorescein (FAM), 
1S5 with hexachlorofluorescein (HEX) and WD0513 with 
cyanine5 (Cy5). The FAM and HEX labelled probes were 
quenched using the black hole quencher 1 (BHQ1) while 
the Cy5 probe was quenched with BHQ3. 

Amplification and detection were completed on a Light- 
Cycler 480II system equipped with a 96- or 384-well block 
(Roche, Germany). The qPCR reaction was performed as 
follows: an initial denaturation step at 95°C for 5 min 
followed by 45 cycles at 95°C for 10 sec, 60°C for 15 sec 
and amplification at 72°C for 1 sec with single fluores- 
cence acquisition. The reaction was prepared in 10 uL 
using 1 uL of DNA, 5 uL of LightCycler 480 Probe 
Master (Roche, Germany) and the volumes of primers 
and probes listed in (Additional file 1: Table SI). 
Fluorescence was acquired simultaneously for the 3 
fluorophores used: FAM, HEX and Cy5. Data were 
analysed with the absolute quantification module using 
the second derivative maximum algorithm of the 
LightCycler480 II. Each PCR plate included known in- 
fected and uninfected control samples. The qPCR 
assay was able to distinguish between the wMel and 
wMelPop-CLA strains with 100% accuracy on hun- 
dreds of analysed larvae and adults. 

Wing morphometries 

Left wings of mosquitoes (unless damaged) were mounted 
on a slide with Hoyers solution [32]. Each wing was 
photographed via 11.25x magnification with a Nikon 
SMZ1500 (Nikon Corporation, Shinagawa-ku, Tokyo, 
Japan) microscope and camera; photos were digitized with 
tPSUtil and tPSDig2 version 2.16 [33]. Fifteen landmarks 
were selected (Figure 1). Wing centroid size was used as a 
proxy for body size [34], computed as the square root of 
the sum of squares of the Euclidean distances between 
landmarks to the centroid. Landmarks were used for 
shape analyses (see below). Measurements of wing length 
were also taken between the alular notch (landmark 11) 
and the furthermost tip of the wing (approximately at 
landmark 4). 

Repeatability of the landmarks was tested through re- 
peat measurements on 200 wings. We ran a one-way 
analysis of variance (ANOVA) for each x- f y- coordinate 
and centroid size treating individuals as a fixed factor. 
The error then estimates within-individual variation (i.e. 
repeatability of measurements) [35]. Repeatability can 
then be computed as the ratio of among individual vari- 
ance to the sum of variance of among and within indi- 
vidual variance. We obtained a repeatability of > 0.99 for 
all landmarks which is adequate for proceeding with an 
analysis of shape [35] . 



Statistics 

All size data were analysed in R3.01, while all shape ana- 
lyses were performed in MorphoJ, but graphed in R. Wing 
centroid size and thorax length data were tested for 
normality using Shapiro- Wilks tests. Variances were 
compared to test for deviations from homoscedasticity. 
Pairwise-differences in means were tested via student t- 
tests, or unequal variance Welch £-test or non-parametric 
Mann- Whitney U tests, depending on whether there was 
unequal variance between groups and/or a breach in nor- 
mality assumptions in one or both groups. ANOVA was 
only performed when data were normally distributed and 
groups had equal variances. CoVs were compared follow- 
ing Miller [36]. For analyses that required multiple com- 
parisons, ^-values were adjusted via the Dunn-Sidak 
procedure. 

Procrustes superimposition was first performed on 
landmark coordinates to standardise wings to one unit 
centroid size, to remove orientation and location effects 
[37,38]. Covariance matrices were generated for each 
superimposed dataset to allow exploration of variation 
via Principal Component Analyses (PC A). To maximise 
separation between groups of interest, we used canonical 
variate analysis (CVA). If within group covariance matri- 
ces were very different between groups or only two 
groups are considered, pairwise group separation was 
performed via discriminant function analysis (DFA). 
Pairwise comparison tests were performed via permuta- 
tion tests on Procrustes and Mahalanobis distances. Pro- 
crustes distance is defined as the square root sum of 
squares of the Euclidean distance between the individual 
Procrustes superimposed shape when compared to the 
superimposed average shape. Mahalanobis distances are 
the square root of the distances squared between the 
superimposed individual to the mean shape that are 
standardised by the covariance matrix of the distance 
variables. In situations where the ^-value differed dras- 
tically due to anisotropy of the variation within groups 
and other factors, we took a conservative approach in 
evaluating the statistics. 

Daily trapping rates across time of collection in DSTs 
were compared via Kruskal-Wallis tests (non-parametric 
ANOVA), with Mann- Whitney U tests used for all pair- 
wise comparisons. 

Results 

Morphometric traits and fitness of released mosquitoes 

Infected mosquitoes caught in the field within 10 days after 
the first release did not differ significantly from semi-field 
cage reared mosquitoes for wing centroid size (CS) or 
thorax length (TL) (see Table 1 for data, CS: t = 0.9838, 
df = 239, p > 0.32; TL: t = 1.4278, df = 263, p > 0.15). Wing 
length (WL) was also not significantly different between 
the two groups (field cage vs field caught: 2.95 vs 2.94 mm, 
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Figure 1 Position and order of landmarks on Aedes aegypti wing. 
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Mann- Whitney p > 0.45). However, wing/ thorax ratio was 
significantly different (see Table 1 for data, Mann- Whitney 
U, p- 0.043); field cage females had a lower wing/thorax 
ratio, i.e. higher wing load than those recollected from the 
field, though they are a subset of the same population. Co- 
efficients of variation for all three morphometric traits 
were similar (Z < 1.9, p > 0.05). 

Since released infected mosquitoes were reared under 
high nutrition, we expected released wMelPop-CLA-in- 
fected female mosquitoes to be larger than uninfected 
female mosquitoes. This is reflected in wing size and 
thorax length (Table 1), including wing length, as com- 
parisons were significant (CS: Welch-£ = 16.36, df = 
381.1, p< 0.001; TL: Welch-/: = 17.10, df= 440.7, p< 
0.001; WL: Mann- Whitney U, p < 0.001). Wing/thorax 
ratio (Table 1) of infected females also differed from 
uninfected females (Mann- Whitney U, p < 0.001), with 
infected females having a higher wing load than unin- 
fected females. The coefficient of variation (CoV) was 
significantly lower (less than half in the case of the size 
traits) in infected females when compared with unin- 
fected females (Z > 5, p < 0.001) as expected since field- 



emerging mosquitoes would have been reared under a 
range of conditions affecting size. 

PCA on shape landmarks for all DST and field cage 
data showed an even distribution across the first four 
principal components (PC) accounting for 62.4% of total 
shape variation (Figure 2). Around 25% (PCI) of vari- 
ation can be explained by landmark 1 changing in the 
opposite direction to landmark 14 and 15; 16.5% (PC2) 
of variation is characterised by landmark 1, 14 and 15 
tending towards the right as the outer landmarks tend 
towards the left; 12.6% (PC3) of variation describes an 
expansion of outer landmarks relative to inner landmark; 
8.1% (PC4) of variation is explained by landmark 12 
moving closer to the other inner landmarks while land- 
marks 7 to 10 move closer inwards. 

A CVA (Figure 3) with DST infection status and field 
cage treated as separate groups indicated some difference 
between infected and uninfected caught in the field based 
on the first canonical variate (CV1). Along CV2, infected 
females caught in the field were differentiated from field 
cage females. These differences were statistically significant 
by permutation tests (10000 replicates) on Mahalanobis 
and Procrustes distances between groups (p < 0.0001). 

Across time there were few changes in the measured 
morphometric traits. Data were separated by infection 
status because of significant heteroscedasticity among 
groups, involving a relatively higher variance in unin- 
fected mosquitoes. Among the uninfected females, there 
was no evidence of differences for all morphometric 
measurements (CS: F 5> 239 < 0.55, p > 0.74; WL: F St 239 < 
0.51, p > 0.77; thorax: F 5) 280 < 0-32, p > 0.9; wing/thorax 
ratio: Mann- Whitney U statistics for pairwise compari- 
sons, p > 0.05). Based on all pairwise comparisons using 
Mann-Whitney U statistics, infected females in DSTs 
also did not exhibit changes in any of the morphometric 
measures (p > 0.05). Pairwise comparison tests of Maha- 
lanobis and Procrustes distances revealed no shape dif- 
ferences with time of collection, except in the case of 
infected females from 9 January and 16 January collec- 
tions (p < 0.001). 



Table 1 Sample size, n, mean/median, standard deviation, SD, and coefficient of variation, CoV of morphometric traits 
by groups of Aedes aegypti from double sticky traps and field cage 



Wing centroid size (mm) Thorax length (mm) Wing size/thorax ratio 





n 


Mean 


SD 


CoV 


n 


Mean 


SD 


CoV 


n 


Median 


SD 


CoV 


Field cage 


91 


3.211 


0.143 


4.46% 


100 


1.468 


0.070 


4.74% 


91 


2.180 


0.047 


2.16% 


Cairns Jan 2012 


























Infected (1) 


92 


3.206 


0.152 


4.75% 


97 


1.465 


0.080 


5.43% 


92 


2.194 


0.054 


2.47% 


Infected (2) 


150 


3.192 


0.148 


4.64% 


165 


1.454 


0.078 


5.35% 


150 


2.199 


0.057 


2.59% 


Negative 


245 


2.827 


0.294 


10.39% 


286 


1.265 


0.157 


12.39% 


244 


2.242 


0.089 


3.99% 



Infected (1) excludes wMelPop-CLA-infected females caught on 16 January because the excluded date would have included individuals from the second release. 
Infected (2) includes all females trapped during the DST trapping period. 
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Figure 2 Shape variation in field uninfected, infected field cage and infected field released A. aegypti. Principal components showing the 
variation in shape for which PC1 explained the greatest variation followed by PC2, PC3, PC4 and so on. (A) PC1 vs. PC2 and (B) PC3 vs. PC4. 
Means and standard deviations of each group and ellipse outlining 90% of data, showing the high degree of overlap in variation for all groups 
within the double sticky trap (DSD and field cage data. Both wMelPop-CLA-infected and uninfected mosquitoes were caught in DSTs. (C) Shape 
variation in PCI, PC2, PC3 and PC4. The procedure over-exaggerates the variation, thus the shape variations were scaled down 10 times. 



Assessing oviposition success of released, infected female 

A total of 165 wMelPop-CLA-infected females and 287 
uninfected females were trapped from 4th January to 
16th January and these mosquitoes provided an oppor- 
tunity to evaluate the oviposition success of infected fe- 
males before females from the ensuing release were 
captured. The rate of capture of uninfected females for 
each house on a per day basis was stable (Kruskal-Wallis 
X 2 = 4.08, df = 5, p > 0.54). On 9 January, there was a 
peak in infected females captured in the DSTs, with 
numbers comparable to those for uninfected females 
(Figure 4). The number and rate of trapped infected fe- 
males dropped substantially before moving up again on 
16th January as a result of the second release on 



11 January. After excluding data for the 4th, 6th and 
16th January (infected mosquitoes that were not ovipo- 
siting or not from the first release), traps caught 96 in- 
fected females vs. 222 uninfected females, or around 
30% of the ovipositing females were infected. 

We estimated possible contributions of the first release 
to the DST collections beyond the 16th January by 
extrapolating from data collected 9th to 13th January. 
Based on a logistic growth model or exponential decay 
function, we estimated that on the 16th January there 
was a contribution of 1 (1%) or 9 (8.3%) wMelPop-CLA- 
infected mosquitoes from the first release, respectively 
(Additional file 2: Figure SI and Additional file 3: Figure 
S2, Additional file 4: Supporting Text). Assuming 
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CV1 vs CV2 



Status 
-I— Field Cage 

wMelPop-infected 
Negative 





3 


rr 










wMelPop-CLA-infected - uninfected 




Field cage infected - released infected caught in field sticky traps 

Figure 3 Shape differentiation between field uninfected, 
infected field cage and infected field released A. aegypti. 

Canonical variate analysis (CVA) treating infection status from DSTs 
and field cage as groups. (A) Canonical variate plot of mean and 
standard deviation, including 90% data ellipse for each group using 
first and second canonical variates (CV1 and CV2) that explains the 
greatest differentiation between groups. CV1 explains 61.9% of 
variation among groups while CV2 explains 29.8%. (B) Shape 
changes from wMelPop-CLA-infected females in DSTs to uninfected 
females in DSTs, and shape changes from field cage females to 
sticky trapped w/MelPop-CLA-infected females. All shape changes 
were magnified 10 times for easier visualization. 



indefinite mosquito survival, the two models predicted 
that 1 (1%) or 15 (13.8%) mosquitoes would have been 
caught in the first gonotrophic cycle after the 13th 
January respectively (Additional file 2: Figure SI and 
Additional file 3: Figure S2, Additional file 4: Supporting 
Text). This means that 86-99% of gravid females in the 
first gonotrophic cycle will seek an oviposition site 
within the first 10 days of release (or before adult fe- 
males reach an age of 13 days), or 52-60% within five 
days after release (before an age of 8 days). 



Effects of infection on size and shape in generations 
subsequent to release 

BGS-trap data from November-December 2012 were sep- 
arated by sex because males and females were significantly 
differentiated for wing size and wing length (males smaller 
than females: median CS 2.19 vs. 2.74 mm; Mann- 
Whitney U, p < 0.001) and shape (pairwise comparison of 
Mahalanobis and Procrustes distance, p< 0.0001). PCA 
yielded a first principal component, PCI, which showed a 
heterogeneous distribution of variation among sexes and 
which can be further differentiated via DFA (Additional 
file 5: Figure S3). The separation in shape is distinct 
enough (more than size) to show that there was an indi- 
vidual wrongly sexed as a female. PCA on female and male 
data indicated a homogeneous distribution of variance 
(not shown), with noticeable variation attributed to the 
separation of landmark 1 from landmark 14 and 15 in 
both PCI and PC2. 

The overall mean wing size of uninfected and naturally- 
reared wMelPop-CLA-infected mosquitoes was not sig- 
nificantly different for both males (CS: Means 2.20 vs. 2.16 
mm, t= 1.22, df = 193, p > 0.22) and females (CS: Means 
2.76 vs. 2.73 mm, * = 0.73, df = 179, p > 0.47). Coeffi- 
cients of variation for wing size also did not differ 
between uninfected and infected mosquitoes for both 
males (CoV: 9.52% vs. 8.31%, Z<1.14, p > 0.25) and 
females (CoV: 10.32% vs. 11.93%, Z < 1.27, p> 0.21). The 
same pattern was evident if wing length was considered. 
There were no differences in overall shape for both 
males and females between uninfected and wMelPop- 
CLA-infected mosquitoes when running permutation 
tests comparing Mahalanobis distance (Male: p > 0.38; 
Female: p > 0.29). 

The amount of size variation measured by the CoV ap- 
pears to fluctuate (See Table 2), but CoV was not signifi- 
cantly different within and among both infected and 
uninfected mosquitoes (separated by sex) (Z<1.96, p> 
0.05). Only the CoV of uninfected and infected females 
from the second collection was significantly different (Z > 
1.96, p < 0.05) (see CoVs in Table 2). 

We note that the average temperature in the 20 days be- 
fore each trapping period was approximately 1.5-2.0°C 
higher in the two December sessions (27.5-29.5°C) than 
the three November sessions (25.8-27.9°C), so we decided 
to compare individuals grouped by month of trapping. In 
pairwise comparisons, uninfected mosquitoes trapped 
in November tended to be larger than those collected in 
December. Males from the first and third trapping session 
(7th and 21st November) were both significantly larger 
than males from 12th December and 19th December 
(p < 0.05 after correcting for multiple comparisons). Simi- 
larly, females from the first and third session (7 and 21 
November) were both significantly larger than females 
from 19th December (p < 0.05 after correcting for multiple 
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Figure 4 Number of females trapped and associated trapping rate in DSTs. (A) Total number of females caught in DSTs over six collections 
(4th, 6th, 9th, 1 1th, 13th and 16th January). (B) Female trapping rate for each house per day for all six collections with 95% confidence intervals 
(number of houses = 50). Both graphs were separated into uninfected and infected individuals. Red dotted vertical lines indicate day when DSTs 
were first deployed while green vertical lines are the first two releases of i/vMelPop-CLA-infected mosquitoes in Machans Beach. 



comparisons), although this was not the case for females 
from 14th November. However, on the 14th November 
the density was much higher; 130 mosquitoes were 
trapped vs 78-97 during the other trapping sessions. Con- 
versely to uninfected mosquitoes, infected mosquitoes did 
not show any clear distinction in wing size and length be- 
tween collections after accounting for multiple compari- 
sons (all p > 0.05). 

When comparing infected and uninfected mosquitoes 
from each collection, uninfected males were larger than 
infected males in the first three collections (in November, 



especially 14th and 21st) but there were no differences in 
the December collections (Table 2). Uninfected females 
were larger than infected females only in the first collec- 
tion, and smaller than infected females in the last collec- 
tion (Table 2). Here, we emphasise that there appears to 
be a trend for decreasing size between November and 
December in uninfected mosquitoes but no significant 
trend in infected mosquitoes. 

Permutation tests on pairwise comparison of Mahalanobis 
and Procrustes distance for shape suggested differences 
between the three November collections and the two 



Table 2 Size measures and statistics of BGS-trap samples of Aedes aegypti separated by collection, sex and infection 
status 



Wing centroid size (mm) 



Uninfected 



wMelPop-CLA-infected 



Sex 


Collection 


n 


Mean 


SD 


CoV 


n 


Mean 


SD 


CoV 


Comparison 


Females 


1 


16 


2.96 


0.27 


9.16% 


12 


2.63 


0.18 


6.91% 


p<0.01 




2 


31 


2.70 


0.24 


9.06% 


17 


2.72 


0.43 


15.96% 


NS 




3 


23 


2.88 


0.27 


9.40% 


8 


2.79 


0.30 


10.66% 


NS 




4 


25 


2.74 


0.28 


10.07% 


10 


2.70 


0.35 


12.95% 


NS 




5 


32 


2.65 


0.27 


10.39% 


7 


2.89 


0.17 


5.86% 


0.05 >p> 0.01 


Males 


1 


34 


2.29 


0.18 


7.82% 


18 


2.19 


0.16 


7.46% 


0.10 >p> 0.05 




2 


35 


2.21 


0.21 


9.42% 


16 


2.06 


0.17 


8.07% 


0.05 >p> 0.01 




3 


13 


2.33 


0.09 


3.77% 


9 


2.19 


0.16 


7.46% 


p<0.01 




4 


30 


2.13 


0.21 


9.67% 


5 


2.29 


0.26 


1 1 .20% 


NS 




5 


29 


2.07 


0.20 


9.82% 


6 


2.15 


0.13 


5.88% 


NS 



Collections 1 to 3 are in November 2012, 4 and 5 are in December 2012. 
n: sample size, SD: standard deviation, CoV: Coefficient of variation. 

Wing size pairwise comparisons were made between uninfected and i/i/MelPop-CLA-infected from the same collection. NS = Non significant, we also reported the 
p-value if less than 0.1. 
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December collections for both uninfected males and fe- 
males (0.05 >p> 0.0001). Infected males and females did 
not differ between collections for the same test. When col- 
lections were grouped into month of collection, both unin- 
fected males and females showed significant difference in 
distances between November and December collections 
(p < 0.001) but not among infected individuals (p > 0.05). 
Landmark 1 was closer to landmark 13 and 15 in Novem- 
ber uninfected mosquito samples (Figure 5). There was no 
clear evidence of differences in distances between unin- 
fected and infected mosquitoes at any time point. 



Discussion 

Released mosquitoes containing the wMelPop-CLA in- 
fection were found to be large with a low variance com- 
pared with uninfected mosquitoes from the field. A 
similar result was obtained with released wMel- infected 
mosquitoes [11]. Within the released mosquitoes, there 
was no evidence for direct selection on size because size 
of field cage reared and captured mosquitoes were simi- 
lar. Again, this is largely consistent with the results for 
wMel [11] and may simply reflect the narrow variance in 
size of the released mosquitoes. Size also did not differ 



Histogram of discriminant values 
(Females) 



Histogram of discriminant values 
B (Males) 




£ CD 




Discriminant values 



Discriminant values 




Females (December - November) 




Males (December - November) 

Figure 5 Shape differentiation in uninfected mosquito caught in November and December. (A) Discriminant value histogram for 
uninfected females comparing between months of BGS-trapping. (B) Discriminant value histogram for uninfected males. (C) Shape change for 
uninfected females from December (points) to November, and shape change for uninfected males from December to November. All shape 
changes were magnified 10 times. 
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among collection time points, suggesting no bias for 
time to successful blood feeding. 

In contrast, there was evidence that the wing/thorax 
ratio differed between the cage and field samples, sug- 
gesting that direct selection had occurred on this trait 
and/or that wing/thorax ratio had influenced ability to 
locate the DSTs. In insects, wing/thorax ratio can con- 
tribute to flight ability and dispersal For instance, field 
releases of Drosophila showed that under field condi- 
tions, flies that have dispersed relatively further have a 
higher wing/ thorax ratio [39]. In the case of released 
mosquitoes, individuals with a higher ratio (lower wing 
load) may achieve relatively greater success at locating 
breeding sites or exhibit a relatively higher survival. Des- 
pite average dispersal of A. aegypti being low and < 200 
m [40,41], the ability to disperse is likely to be an im- 
portant fitness determinant in behaviours such as "skip 
oviposition", host-seeking and oviposition site seeking 
[42-46]. The significant difference in shape may also be 
related to flight ability, but it is unclear if there are phys- 
ical effects of subtle wing shape differences on flight 
ability. 

The DST results suggest that about 30% of females 
caught were wMelPop-CLA-infected. This is substan- 
tially lower than the proportion of the population con- 
sisting of infected mosquitoes, which was estimated to 
be 53.0% (SA Ritchie et al, unpublished) based on the 
approach outlined in [28]. These results suggest 38.3% 
successful oviposition in wMelPop-CLA-infected relative 
to uninfected females (Supporting Text). It thus appears 
that released wMelPop-CLA-infected females perform 
poorly on blood-feeding, female mating success, and/or 
oviposition site seeking in comparison to uninfected 
mosquitoes. Although this Wolbachia infection is known 
to cause severe age related fitness reduction [2,9,10], the 
cohort used in this estimation is relatively young (five to 
13 days old), thus we cannot suggest that this is also 
age-related. 

We estimated that most of the wMelPop-CLA-infected 
mosquitoes (86-99%) would have oviposited before 10 
days after the release. This meant that most mosquitoes 
took a blood meal within 7 days after the release, sug- 
gesting a daily blood feeding success of 46-50%. Daily 
biting rate per mosquito in one study was approximated 
to 0.63-0.76 [47] but this accounted for multiple blood- 
feeds via histologic methods. Our estimates are similar 
to those obtained in previous DST studies on Aedes 
[45,48]. For instance, Marini et al. [48] suggested that 
50-75% of marked-released, unfed Aedes albopictus fe- 
males became gravid in the first five days after release. 

The wMelPop-CLA-infected mosquitoes collected well 
after releases were terminated (i.e. subsequent genera- 
tions breeding in the field) had a similar size, wing/ 
thorax ratio and shape to uninfected mosquitoes. This is 



consistent with the absence of any substantial differences 
in naturally reared wMel- infected mosquitoes when 
compared with uninfected field mosquitoes [11]. How- 
ever, wing size in wMelPop-CLA-infected mosquitoes 
was relatively constant, whereas the wing size of unin- 
fected mosquitoes generally decreased in response to an 
increase in density and/or especially temperature, as ex- 
pected [49-52]. Laboratory experiments suggest that 
wing size of wMelPop-CLA-infected mosquitoes may 
not change much even when mosquitoes develop for 
longer and emerge late under high density conditions 
(PA Ross et al, unpublished). As large size tends to be 
associated with higher fitness [11,53,54], wMelPop-CLA- 
infected mosquitoes may be at a particular disadvantage 
during periods where environmental conditions result in 
a longer larval development time (i.e. lower tempera- 
tures), because the costs of slow development are then 
not countered by a relatively larger size. 

Conclusions 

The findings of this study suggest that wMelPop-CLA- 
infected mosquitoes have low fitness in the field. The 
lower oviposition success rate and smaller size of in- 
fected mosquitoes under some environmental conditions 
mean that the unstable equilibrium that has to be 
exceeded for successful invasion may be high for this in- 
fection even in the wet season [55]. In the dry season, 
additional fitness costs associated with egg quiescence 
are likely to increase the unstable point even further, 
and even if invasion is successful during the release, the 
infection may not persist [3,7] as was observed [16]. The 
lower wing/thorax ratio of released mosquitoes also sug- 
gests that dispersal ability may be limited in released 
mosquitoes. These effects on mosquito fitness and 
morphology were successfully detected using double 
sticky traps (DSTs) when coupled with mosquito moni- 
toring based on BGS -traps. We recommend that future 
biocontrol releases using Aedes aegypti or related species 
should include a deployment of DSTs or the new gravid 
Aedes trap [56] for evaluating mosquito quality. This ap- 
proach provides a rapid assessment of fitness under field 
conditions, which is likely to be more informative in pre- 
dicting field performance than laboratory-based mea- 
sures of fitness. 

We also make a number of recommendations for fu- 
ture releases with this strain. (1) Release programs 
should aim to release a high frequency of wMelPop- 
CLA-infected mosquitoes (comprising >50% of the exist- 
ing population). If high release numbers are not possible, 
successful invasion may depend on suppression of all life 
stages of the existing population before releases are initi- 
ated [57]. This will be particularly important when 
releases are undertaken in areas where estimated 
population sizes are large [58]. (2) Releases should take 
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place over a long period of time. An extended release 
period helps to ensure that Wolbachia frequencies will 
exist over an unstable point for some time. (3) Nutri- 
tional regimes that lead to optimal wing/thorax ratios 
but maintain large adult size need to be explored. This 
will require a detailed understanding of the relation- 
ship between nutritional components and the asym- 
metric change in thorax length and wing size. By 
increasing the overall amount of food available, wing/ 
thorax ratio will probably always decrease since thorax 
length appears to increase at a relatively faster rate 
than wing length with increased food [11]. Because of 
these factors, invasion of wMelPop-CLA is therefore 
only likely to occur in relatively isolated populations 
and/or when releases are tied to another modality, 
such as vector control to reduce the uninfected popu- 
lation, to assist in the spread of the infection [57]. The 
wMelPop-CLA strain may need to be modified through 
multi-generational artificial selection to alter traits 
such as egg desiccation resistance or insecticide resist- 
ance prior to release [13]. 

Additional files 



Additional file 1: Table SI. Primers/probes and optimal concentrations. 

Additional file 2: Figure SI. Logistic growth model on cumulative 
number of wMelPop-CLA-infected females trapped. Cumulative number 
of infected female mosquitoes trapped in double sticky traps (DST) versus 
number of days after the first release. The estimate is based on numbers 
of infected females caught on 9th, 1 1th and 13th January 2012, with the 
first two time points (4th and 6th January) assumed to be close to zero. 
The estimated curve is most likely an underestimate as we lack data for 
any potential lag phase. 

Additional file 3: Figure S2. Exponential decay fitted on wMelPop-CLA- 
infected female mosquito trapping rate. Observed and estimated infected 
female daily trapping rate per house over time, based on rates for 9th, 1 1th 
and 13th January 2012. 50 houses were involved in this study. 

Additional file 4: Supporting Text. Information on the calculation of 
number of expected females that will be caught in the first gonotrophic 
cycle beyond 13th January 2012 (9 days after the first release) based on 
the two models described in Figures S1 and S2. This is followed by the 
derivation of the relative oviposition success rate of released wMelPop- 
CLA-infected females to field uninfected females. 

Additional file 5: Figure S3. Shape differentiation between males and 
females. Discriminant values from discriminant function analysis of wing 
shape of BGS-trap mosquito samples based on sex. One 'female' recorded 
a discriminant value of 18.48, which is an outlier compared to the other 
females and likely to be a misidentified individual. 
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CLA: Mosquito cell-line adapted; DST: Double sticky traps; BGS: BioGents 
sentinel; FAM: 6-Carboxyfluorescein; HEX: Hexachlorofluorescein; 
Cy5: Cyanine5; BHQ: Black hole quencher; CoV: Coefficient of variation; 
PCA: Principal component analysis; CVA: Canonical variate analysis; 
DFA: Discriminant function analysis; CS: Centroid size; TL: Thorax length; 
WL: Wing length. 
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